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The structure and electronic density of states in layered LnFeAsOi-zF^ (Ln=La,Sm; x=0.0, 
0.125, 0.25) are investigated using density functional theory. For the x=Q.Q system we predict a 
complex potential energy surface, formed by close-lying single-well and double-well potentials, which 
gives rise to the tetragonal-to-orthorhombic structural transition, appearance of the magnetic order, 
and an anomaly in the specific heat capacity observed experimentally at temperatures below ~140- 
160 K. We propose a mechanism for these transitions and suggest that these phenomena are generic 
to all compounds containing FeAs layers. For x >0.0 we demonstrate that transition temperatures 
to the superconducting state and their dependence on x correlate well with the calculated magnitude 
of the electronic density of states at the Fermi energy. 

PACS numbers: 74.25.Jb, 61.50.Ah, 75.25.-l-z 



The discovery of a new superconductor LaFeAsOi-^F^ 
with a high transition temperature (T c =26 K) [T] has 
triggered a global search for other Fe-based alternatives 
to Cu-based superconductors, which have dominated the 
field since their discovery in 1986 [2]. Substituting As, 
Fe, and La for other pnicogens |3], transition metals |2] 
and lanthanides [5] , respectively, applying external pres- 
sure [BJ, and optimizing the doping level have pushed the 
T c to 54.5 K [7]. However, since then its value seems to 
have saturated. As doubts have been expressed that T c 
can be raised any further [Hj, it became apparent that 
generic guiding principles for the T c optimization need 
to be developed. 

LaFeAsO is a member of the layered Fe-pnicogens, in 
which FeAs sheets are separated from each other by spac- 
ers such as layers of ionic oxide, e.g. LaO in LaFeAsO, 
[Fig. llTa)] or metal atoms, e.g. Ba in BaFe2As2 [9]. 
In spite of the difference in the nature of the spac- 
ers, Fe As-based materials show surprising similarities 
in the temperature dependence of their structural pa- 
rameters, anomalies in the electric resistance and spe- 
cific heat capacity, and in their magnetic properties 

(e.g. nana in nana). 

A series of theoretical and computational reports ap- 
peared recently describing the electronic properties, mag- 
netic interactions, phonon structure, and the origin of the 
superconductivity in LaFeAsO and related compounds 
(e.g. [H H3 [ini [T71 EHJ El]). The aim of this work is 
twofold: (i) to develop a model for the phase transitions 
observed in FeAs-based materials and (ii) to investigate 
a correlation between electronic density of states at the 
Fermi energy and the experimentally observed values of 
the T c and its dependence on the doping level. 

The calculations were carried out using density func- 




AF2" 




FIG. 1: (Color online) (a) Structure of lxl LnFeAsO (Ln = 
La, Sm) unit cell, (b) Schematics of several spin configura- 
tions within Fe layers shown for a V2 x v2 supercell. The cir- 
cles show the positions of Fe atoms within the Fe layer. Here 
and below, up and down arrows indicate "up" and "down" 
spins, respectively. See text for details. 



tional theory (DFT), the generalized gradient approxima- 
tion functional PW91 [20] and the projected augmented 
waves method [H] implemented in the VASP code [2"2"] . 
The plane-wave basis set cutoff was set to 600 eV. The su- 
percells containing eight (lxl, Fig.[l|, 16 (V^xV^), and 
32 (2x2) atoms and Monkhorts-Pack grids of 252, 132, 
and 36 fc-points, respectively, were used. For the anal- 
ysis of the electronic structure, the charge-density was 
decomposed over atom-centered spherical harmonics. 

In the first part of the paper we consider the relation 
between configurations of the spins associated with Fe 
3d electrons and the lattice structure. Several ordered 
antifcrromagnetic configurations in a V2xv2 supercell 
are shown in Fig. [lib). In AF1, the spins on the neigh- 
boring Fe atoms are antiparallel. In configurations AF2' 



and AF2" spins are parallel along y- and x-axes respec- 
tively; AF2' and AF2" are equivalent in the case of the 
high-temperature tetragonal (T) phase. 




FIG. 2: Potential energy surfaces for the AF1 and AF2 config- 
urations. Dots correspond to calculated energy values. Open 
circles indicate the spin pairs, which are different in Oi and 
O2 configurations. 



After minimization of the total energies with respect 
to both the atomic positions and the lattice parameters, 
the AF1 configuration maintains the T structure (Ta- 
ble m. Configurations AF2' and AF2" relax to two equiv- 
alent orthorhombic (O) structures Oi and 2 , in which 
the Fe 3d spins along the short Fe-Fe bonds are paral- 
lel and those along the long Fe-Fe bond are antiparallel 
[see Fig. [lie)]- The lattice parameters a, b, and c for 
Oi and O2 relate as ai=bi, 61=02, Ci=C2, and ai>6i 
(see also Fig. [2]). The calculated values for the lattice pa- 
rameters for the low-temperature O-phase agree with the 
experimental data to within 0.4 %. (Table p}. The fer- 
romagnetic configuration is considerably less stable than 
antiferromagnetic ones and is not considered here. 

Integration of the AF2" charge-density within LaO 
and FeAs layers shows that the layers are charged: 
(LaO) +,5 (FeAs)- 15 with S = 0.15 |e|. Thus, one can con- 
sider LnFeAsO as a super-ionic compound, in which ionic 
and ion-covalent bonding within the LnO and FeAs lay- 
ers, respectively, is accompanied by the weak ionic bond- 
ing of these layers. The magnetic moments on Fe atoms 
calculated for AF2 are 1.56 fis (Ln=La) and 1.33 ^b 
(Ln=Sm) . These differ significantly from the values sug- 
gested by Mossbauer measurements (~0.35 hb) [2"3"] . 

To find the energy barrier separating the fully relaxed 
AF2' and AF2" configurations, we calculated the total 
energies Eaf2' an d Eaf2" along the path £1 connect- 
ing Oi and O2 (inset in Fig. [2]). Path l\ is parallel to 
the vector n=(l,-l) in the a-b plane. The Eaf2'(Q) and 



Eaf2"(Q), where Q=a-b, are plotted in Fig.[2j For com- 
parison, we also calculated Eafi (Q) along the path £ 2 1 1 n. 

The calculated values of E lt E 2 , and £3 are 0.005 eV, 
0.025 eV, and 0.15 eV, respectively, for LaFeAsO and 
0.006 eV, 0.026 eV and 0.11 eV for SmFeAsO. We note 
that approximate exchange-correlation functionals, such 
as PW91, can underestimate the values of energetic char- 
acteristics by as much as 100%. More reliable values of 
Ei, E%, and £3, as well as those of Fe magnetic mo- 
ments, can be obtained by applying methods, which in- 
clude exact exchange interaction and allow for coupling 
of different many-electron states [2~4"j . 

At Q—0, AF2' and AF2" have the same atomic struc- 
tures and Eaf2'=Eaf2" , yet, their electronic states are 
different. This situation leads to Jahn- Teller (JT) insta- 
bility [25 and formation of a conical intersection at the 
crossover of the potential energy surfaces (PESs) Eaf2' 
and Eaf2"- Correcting for non-adiabatic behavior near 
the intersection, together with taking into account the 
coupling of many-electron states, introduces an effective 
interaction V, which splits the Eaf2' and Eaf2" into 
a higher-energy single-well potential (Eg) and a lower- 
energy double-well potential (Ed) [25] as shown in Fig. [2] 
We can conservatively estimate that < V < E\. 

The lattice dynamics, described by Eg and Ed, has 
three regimes depending on the temperature (T): 

1. For T < Ei — V, the atoms vibrate near the their 
positions defined by one of the orthorhombic energy min- 
ima of Ed (e.g. Oi). In this case the magnetic structure 
is dominated by AF2' configuration (see Fig. [2]). 

2. For Ei - V < T < E x + V, motion of atoms is 
determined by parabolic branches of Ed, although the 
effect of the barrier separating its energy minima can 
not be neglected. The difference between the average 
distribution of short and long Fe-Fe bonds decreases with 
increasing temperature, which corresponds to a gradual 
transition from O to T symmetry. Magnetic order is lost 
because the Fe spins adjust themselves to the momentary 
local atomic structure, so as the spins are parallel for 
Fe atoms forming short Fe-Fe bonds and anti-parallel 
otherwise. In other words, thermal fluctuations of Fe-Fe 
bond lengths cause reorientation of Fe spins (Fig. [2| . 

3. For T > Ei + V, the lattice dynamics is determined 
by parabolic branches of Eg and Ed and the effect of 
the barrier in Ed can be neglected. The lattice has the 
T-symmetry. There is no magnetic order because the 
orientation of the spins changes according to the local 
atomic structure, as described in 2, and also due to cou- 
pling of electronic states of Eg and Ed ■ 

Experimental observations of the structural and mag- 
netic phase transitions in LaFeAsO (e.g. [TOl [HI l2"B"] ) 
suggest that the T— s-O transition takes place gradually, 
with the Q—a~b order parameter exhibiting two kinks at 
T m ax (~160 K) and T min (-140 K), and that the mag- 
netic phase transition occurs at T m i„ or slightly below 
it. In addition, specific heat capacity displays two peaks, 



which also seem to coincide with T max and T m i n [TT| [T3] . 
Similar data have been reported for other FeAs-based 
materials [9] . These results are consistent with the model 
for the three regimes of the lattice dynamics outlined 
above, in which two phase transition temperatures T max 
and T min correspond to E\ + V and E\ — V, respectively. 
We can also speculate, that the decrease in the ampli- 
tude of atomic vibrations during T^O transition [27] 
can contribute to the abrupt drop in the electrical resis- 
tivity observed, for example, in [1.. 
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FIG. 3: Density of states for (La,Sm)FeAsO a! Fi_ 2 :. Letters 
O and T refer to the orthorhombic and tetragonal phases, 
respectively. The Fermi energy is at 0.0 eV. 



We now consider the effect of F-doping on the atomic 
and electronic structures of LnFeAsO. The doping pro- 
vides additional electrons to the FeAs layer so as the 
charge distribution becomes (LnO) +l5+a: (FeAs) _ ' 5 ~ :!: and 
the lattice parameter c decreases due to the increased 
inter-layer ionic bonding (Table ll| . This leads to opening 
up of a narrow gap in the N{e) at ~2.5 eV below the Ef 
(not shown). 

We find that the spin-density distribution in the FeAs 
layers is not independent on the arrangement of the F 
impurities. For example, for x=0.25 (\/2 x \/2 cell), the 
spin-down density is localized on a single Fe atom nearest 
to the F~ impurity, while the remaining three Fe atoms 
share the spin-up density. At the doping level of £=0.125 
(2x2 cell), the effect is more subtle. The lowest energy 
state is similar to that of the undoped LnFeAsO: the lat- 
tice structure corresponds to O-symmetry of the \[2 x \/2 
cell and the spin-arrangement is the same as in AF2, al- 
though the values of fip e are reduced to 1.32 (Ln=La) 
and 0.75 \xb (Ln=Sm). We also found a spin-disordered 
state, which has the T-symmetry and is ~7 (Ln=La) and 
~5 (Ln=Sm) meV per Fe atom higher than the ground 
state. Taking into account the generally random distri- 
bution of the F impurities over O lattice sites in realistic 
samples, we suggest that such spin-disordered state real- 
izes in practice. 



For £=0.125 (2x2 cell) we distinguish two sets of non- 
equivalent Ln and As atoms with different values of their 
z-coordinates. The effect of such structure on the lattice 
phonons and on the charge- and spin-density distribu- 
tions needs to be considered separately. 

Finally, we investigate the correlation between the 
structure and doping level and the electronic density of 
states [iV(e)] calculated for the fully relaxed AF1, AF2, 
and doped LnFeAsOi^^Fj. (Fig. [3]). In all cases the 
N(e) near the Fermi energy (ef) is dominated by the Fe 
3d states and the polarization of spin-up and spin-down 
states is negligible. 

In stoichiometric LnFeAsO, Naf2(e) has a pronounced 
depression near £p, while the Nafi(s) has a narrow deep 
minimum separating a steep rise at e < Ep and a peak 
at e>£f [25]. Projecting Nafi(e) on the d-states shows 
that this peak is dominated by d xz and d yz states. The 
same d xz +d yz peaks near sp are evident for the doped 
LnFeAsO (Fig. [3). 



TABLE I: Structural parameters of LnFeAsO i-jcF^ 
(Ln=La,Sm). In all cases crystallographic cell angles 
a, P, and 7 deviate from 90° by less than 0.0005°. Letters E 
and T refer to experiment and theory (this work) respectively. 



details 



A 6, A 



A 



z(Ln) z(As) 



LaFeAsOi-^F, 
0.0 AF1 T 
0.0 AF2 T 
0.0 300 K E [TO] 
0.0 120 K E [TO] 
0.125 T 

0.25 T 
0.14 120 K g [TO] 
SmFeAsOi-zF^, 



5.6873 
5.7305 
5.7031 
5.6826 
5.6829 

5.6873 
5.6844 



5.6899 
5.6672 
5.7031 
5.7104 
5.6829 

5.6831 
5.6844 



8.6185 
8.6948 
8.74111 
8.71964 
8.5630 

8.4859 
8.6653 



0.1448 
0.1433 
0.1413 
0.1417 
0.1560 
0.1452 
0.1562 
0.1477 



0.6383 
0.6438 
0.6517 
0.6513 
0.6405 
0.6394 
0.6410 
0.6527 



0.0 


AF1 


T 


5.5955 5.5918 8.3435 


0.1406 0.6472 


0.0 


AF2 


T 


5.6232 5.5623 8.4142 


0.1396 0.6515 


0.125 




T 


5.5834 5.5834 8.2884 


0.1523 0.6496 
0.1413 0.6479 


0.25 




T 


5.5888 5.5902 8.2046 


0.1529 0.6493 



According to the standard BCS theory of superconduc- 
tivity, the transition temperature T c is proportional to 
(ui) exp[— 1/XN(ef)], where (u>) is a typical phonon fre- 
quency and A is the electron-phonon coupling constant. 
As shown in Ref. [2S], the limitation of T c <40 K, sug- 
gested by Migdal's theorem for BCS superconductors, is 
not justified and, therefore, much higher values of the T c 
can be achieved by optimizing (w), A, and N(ef). We can 
tentatively suggest that (w) and A do not vary strongly 
for FeAs-based compounds, since the conductivity is con- 
fined to the FeAs layers. Then T c can be considered as a 
function of a single parameter N(ef). 



Thus, we consider the correlation between the behav- 
ior of N(e) for e close to Ef (Fig- |3l and experimentally 
observed properties of LnFeAsOi_ x F x superconductors. 
First, we notice that as x increases and £p shifts across 
the d xz +dyz peak, the value of N(sp) increases as well, 
then reaches its maximum and then decreases. The de- 
tails of the peak structure depend of the value of x but 
its general shape is reminiscent of the experimentally ob- 
served dependence of the T c on x (e.g. [T1H^]V 

Furthermore, the maximum of the d xz +d yz peak 
(x=0.0) in SmFeAsO is higher and further away from £p 
than that in LaFeAsO. This correlates with the observa- 
tions that the optimal T c is higher in SmFeAsOi-^F^ 
(46 K, a;=0.15 [5J) than in LaFeAsOi_ x F x (26 K, 
x=0.05-0.12 pQ) and that it is achieved at larger val- 
ues of x. The slope of N(ef) calculated for x=0.125 is 
negative for Ln=La and positive for Ln=Sm, which indi- 
cates that maximum of N(ef) can be found at x <0.125 
for Ln=La and x > 0.125 for Ln=Sm. This is consistent 
with the optimal values of x found for these compounds 
as ~0.11 (Ln=La) p] and -0.20 (Ln=Sm) [12]. 

Finally, we notice that the d xz +d yz peak in LaFeAsO 
is wider than that in SmFeAsO (this is clearly seen for 
x=0.0 and 0.125), which suggests that T c has a stronger 
dependence on x in SmFeAsO as observed in [T^] . While 
these observations say little about the mechanism of the 
superconductivity in FeAs-based materials, they suggest 
that the highest T c can be achieved in those, which have 
the largest magnitude of the d xz +d yz peak close to Ef- 

To summarize, we investigated the PESs for different 
magnetic states of stoichiometric LnFeAsO (Ln=La,Sm) 
and found that the properties of this system are deter- 
mined by two close- lying PESs: a lower-energy double- 
well potential, where each well corresponds to the or- 
thorhombic symmetry, and a higher-energy single-well 
potential of the tetragonal symmetry. This complex po- 
tential energy surface gives rise to three temperature 
ranges, and, therefore, two transition temperatures, and 
can explain the experimentally observed structural phase 
transition, the appearance of the magnetic order, and the 
anomaly in the temperature dependence of the specific 
heat capacity. 

We noticed a correlation between the calculated pro- 
file of N(e) near ef and experimentally observed depen- 
dence of the T c on the dopant concentration x and on 
the type of Ln atom. This correlation can be used for 
computational prescreening of the promising LnFeAsO 
derivatives as well as for predicting optimal dopant con- 
centrations via relatively inexpensive electronic structure 
calculations. 
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